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Revisiting the parametrization of Equation of State of 
Dark Energy via SNIa Data 
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ABSTRACT 

In this paper, we revisit the parameterizations of the equation of state of dark energy 
and point out that comparing merely the x^ of different fittings may not be optimal for 
choosing the "best" parametrization. Another figure of merit for evaluating different 
parametrizations based on the area of the w{z) — z band is proposed. In light of 
the analysis of some two-parameter parameterizations and models based on available 
SNIa data, the area of w{z) — z band seems to be a good figure of merit, especially 
in the situation that the value of x^j^jn for different parametrizations are very close. 
Therefore, we argue that both the area of the w{z) — z band and Xmin should be 
synthetically considered for choosing a better parametrization of dark energy in the 
future experiments. 

Key words: cosmology - cosmological parameters. 



1 INTRODUCTION 

Current obser v ations , such as those of CMB anisotropy 
(Spcrgcl et al\ 20071'). Supern ovae type la (JRiess et al\ 
12004. i2007l: iDayis e^ a/.l l2007t) and larRe s cale structure 
iTegmark et gZ.I |2004| : lEisenstein et al\ 120071 ') , converge on 
the fact that a spatially homogeneous and gravitationally 
repulsive energy component, referred as dark energy, ac- 
count for about 70% of the energy density of Universe. Some 
heuristic models that roughly describe the observable con- 
sequences of dark energy were suggested in recent years, a 
number of them stemming form fundamental physics and 
other being purely phenomenological. However, the nature 
of dark energy still remain mysterious to physicists and 
astronomers although many possible candidates have been 
proposed. Dark energy present in the equations of cosmolog- 
ical dynamics through its effective energy density and pres- 
sure. The ratio of pressure to energy density (the equation 
of state) is very important in the Friedmann equation re- 
gardless of its physical origin. If dark energy is some kind of 
dynamical fluid and its equation of state would likely not be 
constant, but would vary with redshift z or equivalently with 
cosmic time. The impact of dark energy (whether dynami- 
cal or a constant) on cosmological observations can be ex- 
pressed in term of w{z) = p{z)/ p{z) which is to be measured 
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through either the cosmic expansion history J? (2) (obtained, 
for example, using supernova data) or through large-scale 
structure. Therefore, it is sagacious to study the parameter- 
ization of the equation of state of dark energy empirically 
with as few prior assumptions as possible. 

To reveal the nature of dark energy and narrow down 
the candidate list, a very powerful measure is to map out 
the evolution of the equation of state as redshift changes. 
However, in data fltting, we need to parameterize the equa- 
tion of state w^z) in simple form and then constrain the 
evolution of w^z) in terms of the parameters we introduced 
in our parametrization except the case in which w(z) is 
already such as in the quintessence field (IPadmanabharJ 
2OO7I : ICopeland erall l2006': 'Liu fc Li"2 006|: iHao fc Lil[2003a , 
2004 ). phantom field (Caldwell 20021: bao fc Lil (2003b ; 



Liu fc Lil I2OO3I: iLi fc Hao 20 0j), or Cha plygin gas model 
i Kamenshchik et a;.!.2002 : .Ha'ofc Lill2005h . Unquestionably, 
the way we parameterize the equation of state is bound to 
affect our ability to extract information from the data. There 
are many different parameterizations have been introduced 
based on simpli city and the requirement of regular asymp- 
totic behaviors (|Johri fc Rathll2006l . [2007l : |Johrif[200i ). How- 
ever, will these choice of parameterizations give us maximum 
power to extract information from the data? Some analysis 
existing in literatures compared the different parameteriza- 
tions by looking at their corresponding x^, which are jus- 
tified by the generalized likelihood ratio test in statistics. 
But this measure is no longer fair when the "x^ is small but 
the curvature of the likelihood function is very big, meaning 



2 Dao-Jun Liu, Xin-Zhou Li, Jiangang Hao and Xing-Hua Jin 



that the constraints on the parameters are loose although 
the resulting y^ is small. 

In this paper, we introdu ce another fig ure of 

merit in analogou s to Ref. (|Albrecht et al\ 120061 : 

lAlbrecht fc BernsteinI 120071 ). the area of the 'w{z) — z 
band to evaluate the performance of different parameter- 
izations. The justification of this measure lies in that our 
ultimate goal is to constrain the shape of w(yz) as much 
as we can from the data. In our analysis, we will compare 
the parameterizations with identical number of parameters. 
Note that comparing parameterizations with different num- 
ber of parameters based on Akaike's Information Criterion 
(AIC), Bayesian Information Criterion (BIC) or other crite- 
ria is arbitrary in the sense of the criteria one chooses. In a 
sensible Bayesian method, a model is penalized for having a 
larger number of parameters that gives a reasonable fit (not 
the best fit compared to models with more parameters) 
is awarded with increased evidenc e for that model, (see 
for example, Ref. l|Liddle et aLll2006l )). However, this is not 
what we are concerned and the purpose of this paper is just 
to show what is the best way to parameterize the equation 
of state of dark energy for a variety of prevalent models 
with identical number of parameters. Our results show that 
the widely used parameterization, us^z') — wq -\- wxz j (\ + z) , 
is not the one that can tell us most of the information of 
w(jz) in two-parameter parameterization family based on 
the SNIa data. 

Among the many observations that can help to con- 
strain the shape of w{z), SNIa data provide most sensitive 
and straightforward constraints. Therefore, in this paper, we 
will study the effects of different parameterizations on our 
understanding of the evolution of dark energy based on SNIa 
data. 

The outline of the paper is as follows. In section II, we 
discuss the expansion history of the universe and the obser- 
vational variables from the supernova experiments. In sec- 
tion III, two parameterization families and some prevalent 
models of equation of state of dark energy is introduced. 
Throughout the paper, we only consider two-parameter 
models and parameterizations. In section IV, the method 
and results of the analysis is presented. In the last, we con- 
clude with some remarks on the choice of parameterization 
of dark energy. 



2 THE EXPANSION HISTORY OF THE 
UNIVERSE AND SUPERNOVA 

In the framework of standard cosmological model, assuming 
a spatially flat (fc = 0) Friedmann universe, the equations 
governing the expansion of the universe are 

H'' = Hi {^m{\ + zf + (1 - ^m)S{z)\ , (1) 



and 



<? = 



3i(;(z)(I-nM) + l 



(2) 



where H = a/a is the Hubble parameter, q = —a/aH^ is the 
deceleration parameter, SIm = Pm/pc is the cosmic matter 
density parameter and w{z) is the dark energy equation of 
state, which is defined by 

_ Pde 



The dark energy density parameter pde evolves as pde(z) = 
pDEf{z) and 



f{z) = cxp 



l±^^dz' 



1 + z' 



(4) 



To date, Supernovae Type la (SNIa) provide the most di- 
rect indication of the accelerating expansion of the universe. 
For the distant SNIa, one can directly observe their appar- 
ent magnitude m and redshift z, because the absolute mag- 
nitude M of them are assumed to be constant, i.e., SNIa 
are standard candles. The luminosity distance dL{z) is the 
"meeting point" between the observed m{z) and theoretical 
prediction H{z): 



and 



i{z) = M-|-51ogi 



di(z) = (l + z) 



dL{z) 



Mpc 



-1-25 



cdz 
H{Z) 



(5) 



(6) 



PDE 



(3) 



3 PARAMETERIZATIONS OF DARK ENERGY 

Although H{z) is more directly related to the observable 
luminosity distance and then is easier to measure more ac- 
curately, in order to investigate the evolution of dark energy 
with time and the scale factor, constraints on w{z) is essen- 
tially equivalent to that of H{z) and is also crucial for un- 
derstanding the nature of dark energy ( [Huterer fc StarkmanI 
[2003). Since xu(z) is a continuous function with an infinite 
number of values at a finite redshift range, w{z) must be 
modeled using just a few parameters whose values are deter- 
mined by fitting to observations. A merit of using 'w{z) with 
a particular parametrization is to compare the performance 
of different experiments. Note here that no single parameter- 
ization can represent all possibilities for w(z). A reasonable 
parameterization must be accorded with the demand that 
dark energy is important at late times and insignificant at 
early times. 

There exi st plenty of parameterizations for the eq uation 
of state w{a) (|johri fc RathlbOOd [JOOtI : |johrill2004l ) where 
a — {1 + z)~^ , but most of them are purely phenomenolog- 
ical. Maybe, we should consider some of them in the sense 
that they are generalized from the be havior of physically 
motivated sets of models l|Lindeiil2007l l. 

For single parameter models, e.g. w = costard, no dy- 
namics is embodied and can not parameterize the rate of 
change of w and then high flne-tuning is needed. The phys- 
ical symmetry motivated one-parameter models, such as 
topological defects, are not consistent with the observation 
data. More parameters mean more degrees of freedom for 
adaptability to observations, at the same time more degen- 
eracies in the determination of parameters. For models with 
more than two parameters, they lack predictability and even 
the next generatio n of experiments will not be able to con- 
strain stringently IjLinder fc Hutereill2005l ). Therefore, we 
only consider the two-parameter models in this paper. Of 
course, two-parameter models also have limitations; for ex- 
ample, it is hard to describe rapid variation of w^z) in most 
of these models. 

Various two-parameter parameterization approaches 
have been proposed in the literatures. The simplest 
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way to parameterize the rate of change of w is to 
write the first-order Taylor expan sion. This is the linear- 
redshift param eterization (Linear) iJHuterer fc Turner! I2OOII : 
IWeUer fc Albrecht 2002 1 , which is given by 



W = Wo + W\Z. 



(7) 



This parameterization is not viable as it diverges for 
2 >> f and ther efore incompatible wit h the con- 
straints fro m CMB ijCaldwell fc DoranI [2004) and BBN 
(|johrill2002l). The Upadhye-Isha k-Steinhardt parameteriza- 
tion (UIS) (jUpadhve et aZ.ll2005l ) can avoid above problem, 



wq + w\z, if 2: < 1, 



wo + W\ , 



if 2 > 1. 



(8) 



We here mainly consider the following two commonly 
used two-parameter parameterization families: 
Family I: 



Family II: 



W = Wq + Wi 



W = Wo + Wi 



(1 + 2)- 



(9) 



(10) 



where wo and wi are two undecided parameters, n = 
1, 2, • • • . Both of the parameterization families have the rea- 
sonable asymptotical behavior at high redshifts. The case 
with n = 1 in the above parameterization approaches is 
the same as the most pop ular parameterization introd uced 
by Che vallier and Po larski IjChevallier fc Polarskill200lh and 
Linder (|Lindeiil2003l ) parameterization (CPL). This simple 
parameterization is most useful if dark energy is important 
at late times and insignificant at early times. The one with 
n = 2 in Family II is the Jassal-Bagla-Padmanabhan pa- 
rameterization (JBP), which can model a dark energy com- 
ponent that has the same equation of state at the present 
e poch and at high redshifts, with rapid variation at low 
2 (|Jassal et aZ.ll2005l ). 

Another two-parameter parameterization of dark en- 
ergy equation of state, we consider here, comes from the 
dire ct H(z) p arameterization, first suggested by Sahni et 
al. l|Sahni et al. 2003), 



+ Qiil + z) + {I - Q.M - 0.1 - 0.2 



(11) 



which is corresponding to an effective equation of state of 
dark energy (P2) 



w{z) 



-1 



(l + z)[0^+202{l + z)] 

3 [^^222 + {Oi + 2^2)2 + 1 - Om] ' 



(12) 



On the other hand, there are also two-parameter models 
that have direct physical meani ngs. For example, generalized 
Chip lygin ga s model (000 iJBilic e t al! '200 J: IPev et al\ 
12003 : Chimc nto fc LazkodbOOsI : iLiu fc Li,.2005h , which has 
effective equation of state 



w{z) 



'1 + 



a2{l + zf{ai-l) 
2(22 + 32 -f 3)ai -(1 + 2) 



(13) 



4 METHOD AND RESULTS 

We use the Fisher matrix methods to compute the covari- 
ance matrix for the parameters Wi. If the parameters Wi 
gives the true underlying distribution w, then a chi square 
distribution of data values is in proportion to exp(— x'^/2), 
where x^ is determined by 



-E 



(m» 



(14) 



where ai is error of the distance modulus Hi and f/ is 
theoretical prediction to the data. For SNIa data we use 
here, z^**" = M*''(2; {wt}) = m{z) — M. Using Bayes' theorem 
with uniform prior to the parameter, the likelihood of a pa- 
rameter estimate can be described as a Gaussian with the 
same x^j which is now viewed as a function of parameters 
X^ ~ X^{{''^i})- The distribution of errors in the measured 
parameters is in the li mit of high stat istics proportional to 
(see, for example, Ref. (| Albrecht et aZ.ll2006 )) 



exp 



where the Fisher matrix Fij is defined by 



f -i 



a2 2 

d X 
dwidwi 



(15) 



and < • • • > means average over realizations of the data. The 
covariance matrix of the parameters is simply the inverse of 
the Fisher matrix. 



C^j = {F~') 



(16) 



The error on the equation of state w{z) is given by 
(|Nesseris fc PerivolaropoulosI 120051 ') 



Ej dw 






(17) 



where A'' is the number of the free parameters, a^ is function 
of 2, we define the area of w{z) — z band as 



{z)dz, 



(18) 



where the integral interval {zi,Zh) is take n as (0,2^^3;). 

We make use of the full Gold dataset (JRiess et aZ.ll2004l 'l 
(157 data points, < 2 < Zmax = 1.755) and the com- 
bined Essence, Hubb le, SNLS and nea rby supernovae cat- 
alog as compiled by IjDavis et aZ.1 120071 ). for a total of 192 
supernovae, respectively, assuming a fiat universe with en- 
ergy density in matter Om = 0.3. The value of the Hubble 
constant Ho is marginalized analytically. 

The main results are listed in Table [l] and Tabled Ta- 
ble [T] shows the minima of x^, the best-fit values of the 
free model parameters and their standard deviations in the 
two parametrization families and several prevalent dark en- 
ergy models introduced in section [31 The values of the corre- 
sponding areas of the ^1(2) — 2 bands of different parameter- 
izations and models are also shown in Table [l] All of these 
quantities are worked out by using 157 Gold SNIa data. As 
a comparison, Table[2]shows the results of the same physical 
quantities based on the newly compiled 192 SNIa data. To 
better explain the results, we plot the w{^z) — 2 bands for pa- 
rameterizations of family I, family II and the selected preva- 
lent models in figure [l] figure [2] and figure |3l respectively. 



4 Dao-Jun Liu, Xin-Zhou Li, Jiangang Hao and Xing-Hua Jin 



In figure |4] and figure [5] we show the portraits of Xmin ~ s 
phase of parameterizations of family I and family II and the 
selected prevalent models in the light of the results obtained 
by using 157 Gold SNIa data and newly compiled 192 SNIa 
data, respectively. 

As is shown in figure |3J it is clear that for parameteriza- 
tions in family I, the minima of Xmin Eind s are coincide with 
each other at n = 1, which corresponds to the widely used 
CPL parametrization. However, in family II, the minimum 
of Xmin does uot coiucide with the one of s. The minimum 
of Xmin is located at n = 4, while the minimum of s is lo- 
cated at n = 3. Note that n = 1 parametrization (i.e. CPL) 
have neither the minimum Xmin nor the minimum s. There- 
fore, if we have to choose a parametrization in family II, 
n = 1 parametrization will not be preferred, and the two 
most competitive parametrizations are n = 3 and n = 4. 
But as the difference between two Xmin is much less than the 
difference between two s, we would prefer the n = 3 param- 
eterization. For the prevalent models we investigated here, 
the value of Xmin of GCG model is much greater than those 
of CPL, UIS, Linear and P2 models, so it is not preferred 
in the sense of data fitting, but it may still be interesting 
because of its physical meaning. 

As an improvement and extension of earlier data, the 
newly compiled data set (JDavis et aZ.ll2007l ) provide us more 
sample data. Compare figure [5] with figure |4l we find that 
although there exit minor changes between the results based 
on the 192 newly compiled SNIa data and those based on the 
157 gold data, the main results remain unchanged. Firstly, 
for parametrizations in family I, the best one is still n = 1 
due to its smallest area of w{z) band, albeit the minimum 
value of x^ for n = 4 parametrization is slightly smaller 
than that of n = 1. Secondly, as is also shown in figure |4j 
parametrizations in family II, as a whole, have both smaller 
areas of the w{z) band and lower minimum value of x^ than 
those in family I, and among the parametrizations in family 
II, n = 3 and n = 4 are still the two most competitive ones. 
Finally, among the selected prevalent models, although the 
area of w(z) band of P2 model becomes relative larger, the 
relative locations of these models does not changed signifi- 
cantly, and compared with CPL, UIS, Linear and P2 models, 
GCG model is still not preferred due to its relative larger 
Xmin- The minor difference between figure|3]and figure[S]may 
arise from the data calibration of different data sets among 
which we shall leave in a future work about a comprehensive 
comparison. 



5 CONCLUSIONS AND DISCUSSIONS 

Traditionally, forming the so-called Bayes factor (likelihood 
ratio for frequentists) Bij = L{Mi)/L{Mj), where L{Mi) 
is called likelihood for the model Mi to obtain the data 
if the model is true, is used in co mparison of the cosmo- 
loKical(and/or dark energy) models (|John fc Narlikaij|2002l : 
iNesseris et aZ.ll2005l ). Generally, L{Mi) is dependent on the 
prior probability and the likelihood, which is determined by 
X^, for the model parameters. And when one has no prior to 
the model parameter, everything is determined by x^; which 
is a measure of the fit to the data. However, according to 
the above analysis based on the latest SNIa data, there exist 
lots of cosmological models and parametrizations of dark en- 
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ergy which lead to very similar Xmin- This is especially true 
for the results we obtained based on the newly compiled 
192 data. The best fitting of any parametrization or model 
we consider here is within the la bound of those of others, 
see table [5] or figure [S] Under this circumstance, how do we 
compare them? Or what parametrization approach should 
be used to probe the nature of dark energy in the future ex- 
periments? The Bayes approach only works in the condition 
that fittings of models are distinctly different. When the dif- 
ference of x^ is very small, one should pursue other figures 
of merit. The above introduced area of 'w{z) — 2 band, we 
think, is such a figure of merit and our point is that both x^ 
and the area of w{z) — z band should be synthetically con- 
sidered for choosing a better parametrization of dark energy 
in the future experiments. 

Bearing this point in mind and according to the re- 
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Table 1. The m ininia of x^ and areas of the w{z) band for different models using 157 Gold SNIa data 
JRiess et al\\2004) . 



Model 




wo{ai or fii) 


wi(a2 or Q2) 


s 


UIS 


174.365 


-1.40802±0.255438 


1.70941±0.928001 


1.69918 


Linear 


174.365 


-1.39978±0.249302 


1.66605±0.892594 


2.13868 


P2 


174.207 


-4.16234±2. 62176 


1.67458±1. 06813 


0.895176 


CPL 


173.928 


-1.57705±0.326346 


3.29426±1.69727 


1.62803 


GCG 


177.063 


0.999827±0.00663069 


83.4676±3209.18 


1.58943 


Family I 










n=l 


173.928 


-1.57705±0.326346 


3.29426±1.69727 


1.62803 


n=2 


174.606 


-1.27011±0.192597 


6.20395±3.44687 


2.13084 


n=3 


175.09 


-1.17171±0.154602 


12.8437±7.65729 


2.71702 


n=4 


175.444 


-1.12497±0.138595 


26.4417±16.8124 


3.37756 


Family-II 










n=l 


173.928 


-1.57705±0. 326346 


3.29426±1.69727 


1.62803 


n=2 


173.409 


-1.87262±0.456452 


6.62831±3.29276 


1.14253 


n=3 


172.824 


-2.39635±0.69199 


13.7569±6.65922 


0.731673 


n=4 


172.454 


-3.2745±1. 11672 


28.3698±13.8425 


1.30027 



Table 2. The minima of x'^ and areas of the w(z) band for different models using 192 SNIa data IIDavis et al\ 
I2OO7I: IWood-Vasev ef"aIl l2bo7: Ricss et al. 2007). 



Model 




W(,(ai or r2i) 


w\{a2 or Q.2) 


s 


UIS 


195.412 


-1.1192±0.27732 


0.0485532±1. 17151 


2.21778 


Linear 


195.409 


-1.12628±0.281052 


0.0811196±1. 18901 


2.92336 


P2 


195.382 


-0.591863±1. 71977 


0.166512±0.679303 


3.07726 


CPL 


195.411 


-1.12456±0.331918 


0.0961458±1. 89159 


1.88532 


GCG 


195.529 


1.00055±0.00898214 


95.6168±1553.63 


1.62394 


Family I 










n=l 


195.411 


-1.12456±0.331918 


0.0961458±1. 89159 


1.88532 


n=2 


195.413 


-1.11369±0.21235 


0.135963dz4.91012 


3.1127 


n=3 


195.402 


-1.12407±0.181828 


1.52726±14.2855 


5.13255 


n=4 


195.314 


-1.15242±0.164779 


13.4211±34.8666 


7.01332 


Family-II 










n=l 


195.411 


-1.12456±0.331918 


0.0961458±1. 89159 


1.88532 


n=2 


195.409 


-1.13475±0.412811 


0.203332±3.09753 


1.13799 


n=3 


195.399 


-1.17258±0.546306 


0.631377±5.28164 


0.641504 


n=4 


195.356 


-1.29367±0.787516 


2.28402±9.61853 


0.960275 



2r 
1 

5 
-1 



4 
S 2 


-2 - 









4 

















1 














2 








-0.5 

















/^ 


^^ 




s 




^ 


^ 




S -1 

-1.5 








S -1 

-2 


























-2 








-2 








-3 








0.5 


1 


1.5 


0.5 


1 


1.5 


0.5 


1 


1.5 


0.5 


1 


1.5 




z 








z 








z 








z 






n=3 




6 




n=4 




5 




n=3 




10 




n=4 


































4 
S 2 









% 








5 
S 














^^^ 


^^ 







^ 


^ 




-5 
















~ 










-2 








-10 








-5 









0.5 1 1.5 



0.5 1 1.5 



0.5 1 1.5 



0.5 1 1.5 



Figure 1. w(z) — z band for family I parametrizations. T he left four panels are obtained by using 157 Gold data I Riess et a^.ll2004h and 
the right four are obtained by using latest 192 SNIa data llDavis et ai.ll2007l : IWood-Vasev et aLir2007l : iRiess et a;.ll2007t) . 
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suits presented in the above section, we find that the widely 
used CPL parametrization, which has a very simple inter- 
pretation in terms of the scale factor, is not statistically 
"special" among the two-parameter parametrization fami- 
lies, instead n = 3 in family II, which looks like a variation 
on the CPL parametrization, is more preferred. Note that 
CPL parametrization corresponds to the n = 1 case in both 
family I and family II and if we also take n as a free pa- 
rameter, family I and family II are just two three-parameter 
parametrizations. However, in this work we only consider 
two-parameter parametrizations and n is not treated as a 
free parameter, so n = 1 and n = 3 actually denote two 
distinct parametrizations, in which the evolution of wlz) is 
qualitatively different. 

There is an interesting question that whether the dif- 
ferences among the area of w{z) band are significant enough 
to single out one parametrization. In our opinion, the an- 
swer is somewhat depended on the observational data. As 
far as the data we used here, the differences among the ar- 



eas of w[z) band are so significant that we can pick out 
n = 3 of the family II as the best parametrization among 
the models we consider in this work. However, this does not 
mean that the other parametrizations are completely ruled 
out. For example, the simple CPL parametrization and P2 
model still do well to a certain extent. It should be also 
pointed out that, the differences among the areas of 10(2:) 



band are much more significant than those among Xn 



for 



both 157 gold data and latest 192 data. This fact indicates 
that the area of Vli{z) band is likely to be a good figure of 
merit, especially in the situation that the value of Xmin for 
different parametrizations are very close. 

Generally speaking, the motivation from a physical 
point of view should be at the top priority when we choose 
cosmological models. However, it is perfectly clear that in 
the absence of any compelling dark energy model, the sug- 
gested parametrizations are phenomenological. Then the 
reason why people might prefer a given parametrization is 
because of its simplicity and also because they feel that it 
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allows us to extract useful infomiation for a very large class 
of models, and hopefully the "true" model is one of them. 
Anyway, to estimate the effects of dark energy one needs 
to quantify them and parametrization of w has turned out 
to be an efficient tool in this respect. Therefore, there is a 
subtle balance between motivation from a physical point of 
view and fitting results. To help making decisions in this 
situation, we need to know what is the best achievable fit- 
ting result from various models or parametrizations with 
the same number of parameters. This will serve as a fiducial 
criteria for us to choose a best model. The figure of merit 
introduced in this paper is to help to define what is the best. 
As is well known, besides SNIa observations, there exist 
lots of other experiments probing different aspects of dark 
energy and we will have ma ny more data from these experi- 
ments iJAlbrecht et aZ.l2006l l. However, in terms of constrain- 
ing the evolution of w{z), SNIa approach is the most sensi- 
tive and direct one. Other methods, such as CMB and cluster 
counts, are primarily good for the energy density constraint. 
But it will be advantageous to test all the parametrization 
with all the combined data sets in the future. The current 
analysis in this paper could be directly generalized to the 
case with multi-experiments by maximizing the product of 
the likelihood of each experiment. It is worth noting that the 
best parametrization of dark energy models for SNIa data 
may not necessarily be the best one for other observational 
data. We will report that in a preparing work. 
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